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Abstract 



Anisotropics in the temperature of the cosmic microwave background have 
been detected on a range of scales by several different experiments. These 
anisotropies reflect the primordial spectrum of metric perturbations in the 
early universe. In principle, the largest barrier to a clean interpretation of the 
experimental results is contamination by foreground sources. We address this 
issue by projecting out likely sources of foreground contamination from seven 
separate small-angle and medium-angle experiments. We then calculate likeli- 
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hood functions for models with adiabatic perturbations, first for the amplitude 
of the spectrum while constraining the spectral index to be b = 1, and then 
jointly for the amplitude and spectrum of the fluctuations. All of the experi- 
ments are so far consistent with the simplest inflationary models; for n = 1 the 
experiments' combined best-fit quadrupole amplitude is Q rms _ ps = 18^f fiK, 
in excellent agreement with the COBE two-year data. In (Q rms -ps, n ) space, 
the allowed region incorporating intermediate and small-scale experiments is 
substantially more constrained than from COBE alone. We briefly discuss 
the expected improvement in the data in the near future and corresponding 

constraints on cosmological models. 
98.70.Vc 
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Typeset using REVTgX 



More than any other cosmological observation, measurements of temperature anisotropics 
in the Cosmic Microwave Background (CMB) are a direct probe of the primordial spectrum 
of metric perturbations. Precise measurements will give both the amplitude and spectral in- 
dex of the perturbation spectrum, may allow disentangling scalar from tensor perturbations, 
and will provide information on the origin of the perturbations (whether from inflation or 
some type of topological defect, for example). 

Measurements of CMB temperature anisotropics are notoriously difficult and have only 
recently attained the necessary signal-to-noise ratio for meaningful results. Since the COBE 
satellite made the first anisotropy detection in 1992 [1], nearly a dozen other experiments 
have announced positive detections on a wide range of angular scales at amplitudes of a few 
parts in 10 5 of the background temperature [2]. At such small temperature differences, the 
main experimental obstacle is instrumental noise; technological advances and experimental 
ingenuity have pushed noise levels down by a factor of 10 4 in the past 30 years, and this trend 
will continue for the coming few rounds of experiments. The largest hurdle in principle is 
disentangling the foreground contribution: for a given measurement, how much of the signal 
comes from the blackbody CMB and how much from other sources of microwave radiation? 

Two different techniques are useful for sorting out the foreground. The first is to extrap- 
olate sky maps at other frequencies (e.g. radio maps) to estimate the microwave emission 
in various parts of the sky, and then subtract this from the measured signal to obtain the 
cosmic signal. This process depends on detailed modeling of various sources and involves 
uncontrolled extrapolations over large frequency ranges. In this paper we focus on the com- 
plementary method of using measurements at multiple frequency channels to subtract out 
the non-blackbody piece of the measured signal. Clearly, as measurements improve, a com- 
bination of the two methods will give the most reliable interpretation of experiments; here 
we show how much various measurements may be affected by foreground contamination. 

We analyze the current small-scale and medium-scale anisotropy measurements which 
employ multiple frequencies and for which data is publicly available: MAX [3], MS AM 
[4], South Pole [5], and Saskatoon [6]. The first two are balloon-born packages, while the 
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other two are based on the ground. All measure in either three or four frequency channels. 
Saskatoon measures at the largest angular scale, roughly two degrees, and South Pole is 
slightly smaller; we call these two experiments "medium-angle." MSAM and MAX measure 
at smaller scales, roughly 20 to 30 arc-minute scales; these are referred to as "small-angle." 
The angular scale is more precisely characterized by the window function Wi of each exper- 
iment, defined as follows: If the sky temperature with mean T = 2.735° K is decomposed 
into spherical harmonics as T = T (l + J2i m a imYim) } then the experiment will measure a 
mean temperature variance (in the absence of noise) given by 

\T J ^ 47r m tl/ 

For the purposes of the likelihood analysis presented below, the window function is extended 
to covariances between temperatures measured at different points on the sky. Each experi- 
ment has a unique window function, which incorporates the beam profile, chopping strategy, 
scan pattern on the sky, and other experimental details. We have calculated window func- 
tions for the experiments considered here to an accuracy of around 1%. 

We employ the well-known method of Bayesian likelihood analysis [7] to the data from 
each of these experiments. Likelihood analysis asks how likely is a given set of measurements 
if a particular theory is true. The likelihood function quantifying this probability is 



(2 7 r)-^ 2 
L = — , exp 
VdetC 



--D T C~ 1 D 
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Here D is a data vector of length N; in this case, D contains values of the temperature 
anisotropy measured at N p patches on the sky with N c different channels (either at different 
frequencies or polarization states), such that N p N c = N. The correlation matrix C is the 
expectation value (DD T ) which has two separate pieces, the instrumental noise and the 
theoretical signal. The latter depends on the parameters in the theory being tested. Given 
a likelihood function depending on a two-parameter (p\ } p2) theory, we obtain an allowed la 
region in (^1,^2) by the condition 

J dpi dp 2 C(p\ } p2) = 0.68; £(T) constant, (3) 



where the boundary of the allowed region is denoted by T. For a single parameter, the region 
reduces to an interval and the boundary is its endpoints. 

The above analysis is standard practice. Before applying this analysis to the data, we use 
the multiple frequency channels available to discriminate against foreground contamination 
[8]. In particular, for each experiment we choose one source of foreground judged most 
likely to be a contaminant: for MAX and MS AM dust emission, while for South Pole and 
Saskatoon synchrotron emission. These choices are based on the frequency ranges of the 
experiments: the high frequency experiments are more likely to be sensitive to dust while 
the low frequency ones more sensitive to free-free and syncrotron emission. This component 
is then assigned a given spectral index. For example, the signal due to free-free emission is 
assumed to scale with frequency as S(v) = S(v\)(v\l v) 2 ' 1 . If a given experiment has only 
two frequency channels, then the linear combination 

D = D(v x ) - {v 2 jv x f x D{v 2 ) (4) 

is completely independent of free-free contamination. Note that if the frequencies are closely 
spaced, then D approaches the difference between the data in the two channels. This 
difference is zero for a cosmic signal (expected to be frequency independent), so the signal 
to noise ratio becomes very small in the limit of closely spaced frequencies. All other factors 
being equal, experiments that cover a large range of frequencies are best able to distinguish 
cosmic signal from foreground contaminants. This advantage shows up noticably in our 
analysis. The extension of this method to more channels is straightforward: we always 
choose the A^hannei — 1 linear combinations of the data that are independent of foreground 
contamination from the assumed source. Conversely, an A^-channel measurement can in 
principle be used to project out N — 1 different foreground sources, but in practice the signal- 
to-noise ratios are small enough that substantial redundancy is necessary for reasonably 
significant results. 

The likelihood and foreground analysis described above is completely general. We per- 
form the Bayesian likelihood analysis on the class of theories, based on inflation, with a 
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primordial gaussian spectrum of density perturbations 

((Sp/pf) oc k n 

with k the wavenumber of the perturbation. First we perform the analysis with the spectral 
index fixed at n = 1 and the amplitude Qrms-ps [9] as the only free parameter of the theory, 
and then we consider a two-parameter theory depending on both n and Q rms _ ps . 

This two parameter inflationary theory is well-developed and detailed predictions of 
CMB temperature anisotropics are readily available [10]. The actual CMB temperature 
anisotropics depend also on a variety of cosmological parameters: the Hubble constant 
H = 100/i km sec -1 Mpc -1 , the baryon mass fraction the nature of the dark matter, 
the equivalent mass fraction in a cosmological constant Oa, the tensor perturbation spec- 
trum, and the redshift of reionization zr. All of these parameters may vary in inflationary 
models [11]. Nonetheless, recent work [12] has demonstrated that to within 10% the CMB 
anisotropics depend only on an effective spectral index defined by 

n = n- 0.281n(1.56 - 0.56n) - 0.000364 /2 + 0.26 - 2h^l - fi A ) . (5) 

This relationship assumes the preferred nucleosynthesis relation ft^h 2 = 0.0125 [13] and the 
tensor spectrum conditions nj = n — 1 and r = /C2 ~ 7(1 — n) given by the simplest 
inflation models [14]. The effects of relaxing these conditions has been explored elsewhere 
[15]; it appears that even if these conditions are relaxed the theory can still be parameterized 
in terms of its amplitude and a generalized form of Eq. (5). Thus shifting the value of 
these cosmological parameters only moves to a different place in (n, Q r ms- P s) space. This 
"cosmic confusion," as it has been dubbed [12], means that the CMB will never be a strong 
positive test for any particular set of cosmological parameters. This pessimistic outlook 
was the conclusion of Ref. [12] and has since been often reiterated. On the other hand, we 
want to emphasize that this same effect makes the CMB an extremely powerful negative 
test: if various experiments end up being inconsistent for a particular set of parameters, 
Eq. (5) shows that they rule out all other inflationary scenarios which involve only shifting 
cosmological parameters. 
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Proceeding with the experiment analysis, we first fix h = 1 and consider the conditional 
likelihood for the perturbation amplitude given by the small and medium-angle experiments. 
The best fit value and la range for Qrms-ps are plotted for each of seven data sets in Fig. 1. 
The x-axis shows the region in £-space to which each experiment is most sensitive; the 
angular scale in degrees is given roughly by 100 /£. COBE (which we have not analyzed) is 
sensitive to large angular scales and is displayed for reference at a small value of £ [16]. 

Figure la plots the allowed range of Qrms-ps neglecting the effect of foreground; Fig. lb 
projects out foreground contamination as described above. For MSAM and MAX, projecting 
out foreground makes little difference in the error bars because these experiments cover a 
large range of frequencies [17]. For Saskatoon and South Pole the error bars on Qrms-ps 
become much larger when foreground is projected out [8]; note especially the South Pole 
experiments which cover the narrow frequency range 25 — 35 GHz. This is because the 
effective signal to noise ratio is small after projecting out background from a signal spanning 
a narrow frequency range, as noted above. In every case, the \ 2 P er degree of freedom for 
the most likely amplitude is reduced by projecting out foreground. This suggests that fitting 
just one blackbody component to the measurements is not sufficient: projecting out at least 
one foreground component is essential. 

The agreement between large and small scale observations for the h = 1 model is re- 
markably good. The combined best fit from the seven medium/small scale experiments is 
Qrms-ps = 18lf //K, in complete agreement with the COBE two-year values of 20.4 ± 1.7//K 
[16] and 17.6 ± 1.5//K [18] . The uncertainty on the COBE measurement is not going to 
get much smaller, being dominated by cosmic variance. Cosmic variance is not yet a major 
factor for the smaller scale experiments; the uncertainty on Qrms-ps from these experiments 
is essentially limited only by how well the foreground contamination can be eliminated. 

For a two-parameter fit, we now allow h to vary and find the allowed region in (Q r ms- P s, 
h) space for the small and medium-angle experiments. At large angles, the COBE team 
has performed a similar analysis; we quote their results here [18]. Figure 2 shows the 
regions allowed at a la level by COBE, the medium angle experiments and the small angle 
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experiments after projecting out one foreground component. A large region of consistency 
is currently allowed. Clearly at this stage of the observations, the medium and small-angle 
experiments have uncertainties which are too large to make this a powerful test. The medium 
angle experiments in particular have large uncertainties due to the lack of frequency coverage. 
In the near future the situation will improve [2]; in fact, both Saskatoon and South Pole 
have recently been redone with larger frequency coverage, which should substantially reduce 
the error bars. If the allowed regions for the three classes of experiments eventually do not 
overlap, this will be strong evidence against the inflationary scenario. 

Another noticeable feature of Fig. 2 is the different slopes of the allowed regions for 
the three different types of experiments. The large angle results are least sloped since 
COBE is sensitive to the lowest order multipoles, those closest to the normalization point 
at i = 2. A tilt in the spectrum thus has a relatively small effect on Q rms _ ps . The small 
angle experiments are most affected by a spectral tilt, with their best fit Qrms-ps being 
significantly reduced (increased) for n large (small). The variations in slope are essential 
characteristics of experiments on differing angular scales, ultimately leading to a powerful 
test of inflationary models. 

To conclude, we have analyzed small and medium scale anisotropy experiments by pro- 
jecting out one foreground component. For experiments with wide frequency coverage, this 
procedure does not substantially increase the error bars on the parameters in a theory. All 
current medium and small-angle experiments with multiple frequency channels are consistent 
with simple inflation models normalized to COBE, with error bars comparable to COBE. 
The inflation model predictions are degenerate in the cosmological parameters, making the 
theory of inflation testable since only two free parameters determine the anisotropics. The 
best way to test this theory at present is through CMB experiments at different angular 
scales. Current medium angle data lack the frequency coverage to discriminate effectively 
against foreground, so the test is not yet very powerful. However, in the next three to 
five years, we expect improving signal-to-noise ratios and wider frequency coverage will test 
inflationary models in (Q r ms- P s, n) space at the 10% level. 

8 



ACKNOWLEDGMENTS 



We thank Stephan Meyer and Lyman Page for patient explanations of their experiments. 
This work was supported in part by the DOE (at Chicago and Fermilab) and by NASA 
through grant No. NAGW-2381 (at Fermilab). AK was supported in part by the NASA 
Graduate Student Researchers Program. 



9 



REFERENCES 



[1] G.F. Smoot et al, Ap. J., 396, LI (1992). 

[2] For a recent review of some current experimental results and associated difficulties see 
M. White, D. Scott, and J. Silk, Ann. Rev. Astro. Astrophys. (in press); L. Page, in 
Proceedings of Workshop on Cosmic Background Radiation Two Years After COBE, 
ed. L. Krauss (World Scientific, in press 1994). 

[3] This data is from two scans of MAX3: J. Gundersen et al, Ap. J. 413, LI (1993); P. 
Meinhold et al, Ap. J. 409, LI (1993). 

[4] E. S. Cheng et al, Ap. J. 422, L37 (1994) used two independent chopping schemes 
(MSAM2 and MSAM3) to look at the same sky at two different angular scales. Their 
analysis (presented in the "raw" graph in Figure 1) includes the effects of dust in a 
slightly different way than we do. We utilize the entire MSAM data set; that is, we 
do not drop the regions of data which possibly contain unresolved sources. Hopefully 
further study of the same region will soon resolve whether sources are actually present. 

[5] We analyze two scans taken at the South Pole in 1991: in T. Gaier et al., Ap. J. 398, 
LI (1992); J. Schuster et al, Ap. J. 412, L47 (1993). 

[6] Data taken in Saskatoon in 1993: E. J. Wollack et al, Ap. J. 419, L49 (1993). 

[7] A. C. S. Readhead et al, Ap. J. 346, 566 (1989). 

[8] S. Dodelson and A. Stebbins, Ap. J. 433, 440 (1994). 

[9] In a given model, Qrms-ps can be related to the amplitude of the initial power spectrum, 
A = P(k)/k n . For example, if anisotropics are generated only by scalar perturbations 
to the metric, then 



A=(Q 



rms— ps 




10 



See P.J.E. Peebles, The Large Scale Structure of the Universe (Princeton: Princeton 
University Press, 1980). 

[10] P. J. E. Peebles and J. T. Yu, Ap. J. 162, 815 (1970); M. L. Wilson and J. Silk, Ap. J. 
243, 14 (1981); J. R. Bond and G. Efstathiou, Ap. J. 285, L45 (1984); N. Vittorio and 
J. Silk, Ap. J. 285, 39 (1984); N. Sugiyama and N. Gouda, Prog. Theor. Phys. 88, 803 
(1992); S. Dodelson and J.M. Jubas, Phys. Rev. Lett. 70, 2224 (1993); R. Crittenden, 
et al, Phys. Rev. Lett. 71, 324 (1993). 

[11] We assume = 1 for this analysis, as this is a generic prediction of inflation. Open 
universes have been analyzed by M. L. Wilson and J. Silk, Ref. [10]; L.F. Abbot and 
R.K. Schaefer, Ap. J. 308, 546 (1986); M. Kamionkowski and D.N. Spergel, submitted 
to Ap. J. (1994); M. Kamionkowski, D.N. Spergel, and N. Sugiyama, submitted to 
Ap. J. Letters (1994). 

[12] J.R. Bond et a/., Phys. Rev. Lett. 72, 13 (1994) 

[13] T.P. Walker et a/., Ap. J. 376, 51 (1991). 

[14] R.L. Davis et al, Phys. Rev. Lett. 69, 1856 (1992). 

[15] S. Dodelson, L. Knox, and E.W. Kolb, Phys. Rev. Lett. 72, 3444 (1994) studied the effect 
of relaxing the conditions on the tensor perturbation spectrum, allowing the amplitude 
of the tensor modes to be another free parameter. Previous work varying 0^ and h was 
carried out by J.R. Bond, G. Efstathiou, P. M. Lubin, and P.R. Meinhold, Phys. Rev. 
Lett. 66, 2179 (1991); S. Dodelson and J.M. Jubas, Ref. [10]; K. Gorski, R. Stompor, 
and R. Juszkiewicz, Ap. J. 410, LI (1993). 

[16] K.M. Gorski et a/., Ap. J. in press (1994). 

[17] In Fig. lb, points from both the MAX GUM scan and the //PEG scan are plotted. The 
//PEG scan is apparently highly contaminated by dust emission and is off the scale of 
Fig. la displaying the uncorrected data. Subsequent MAX scans in different regions of 



11 



the sky have confirmed the result from the GUM scan. 
[18] C. L. Bennett et al, Ap. J., in press (1994). 



12 



FIGURES 

FIG. 1. The results of seven medium and small-angular scale CMB experiments, plotted as 
the amplitude Q rms _ ps versus angular scale in multipole number. The effective spectral index 
is assumed to be h = 1. The error bars represent la deviations. Figure la shows the raw 
measurements; Fig. lb shows the same measurements after a likely foreground contaminant has 
been projected out. 

FIG. 2. Likelihood contours plotted in (Q r ms-ps 5 n) space. The overlap region between small, 
medium, and large-angle experiments is currently allowed at a 68% confidence level. 
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